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Abstract 

We present results of conductance-noise experiments on disordered films of crystalline indium ox- 
ide with lateral dimensions 2/xm to 1mm. The power-spectrum of the noise has the usual 1// form, 
and its magnitude increases with inverse sample-volume down to sample size of 2/im, a behavior 
consistent with un-correlated fluctuators. A colored second spectrum is only occasionally encoun- 
tered (in samples smaller than 40/xm), and the lack of systematic dependence of non-Gaussianity 
on sample parameters persisted down to the smallest samples studied (2/j.m). Moreover, it turns 
out that the degree of non-Gaussianity exhibits a non-trivial dependence on the bias V used in the 
measurements; it initially increases with V then, when the bias is deeper into the non-linear trans- 
port regime it decreases with V . We describe a model that reproduces the main observed features 
and argue that such a behavior arises from a non-linear effect inherent to electronic transport in a 
hopping system and should be observed whether or not the system is glassy. 

PACS numbers: 72.20.Ee 72.20.Ht 72.70.+m 
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I. INTRODUCTION 



Conduction noise is a property of essentially all electronic systems. A common form 
of this noise has a 1/ f v power-spectrum with rj of order unity I]. Such a spectrum may 
result from the superposition of many ffuctuators with individual frequencies u>i extending 
uniformly over the observed range [if. When these ffuctuators are independent, the central 
limit theorem mandates that the power-spectrum be Gaussian, and this is presumably the 
generic result in the thermodynamic limit. Therefore when a system with power-spectrum 
of the 1/f type and a non-Gaussian nature is encountered, it seems reasonable to conclude 
that the associated ffuctuators are correlated. 

Identifying the underlying source of the correlations however, is not a trivial task as 
there could be several reasons for non-Gaussian noise, including artifacts, and correlations 
in the noise may appear even when the ffuctuators are basically independent. For example, 
Seidler et al using their dynamical-current-redistribution (DCR) model, argued that a 
non-Gaussian noise may result from features that are peculiar to the process of electrical 
conductivity in a network; Their considerations are especially relevant for inhomogeneous 
systems, and when the fluctuations are strong. Both ingredients are inherent to transport 
in disordered conductors where non-Gaussianity were frequently observed jjj. The archtyp- 
ical system of this class is a hopping system where transport is confined to a percolation 
network. Such a system may be viewed as a random-resistor-network in which each resistor 
Tij is associated with pair of sites i,j that are connected by a hopping process jj]. The wide 
distribution of the r^-'s in the random- resistor- network leads to several unique features of 
electronic transport in such a medium. The most familiar feature is that the current in the 
system is earned by a p—n-netwo* fl Q Q involving a r e,ative, y Sm all nu mb e r of 
'critical' resistors. This percolation problem differs from the classical, 'geometric' scenario 
in two essential accounts; First, the current-carrying network (CCN) of the hopping system 
is temperature dependent - it becomes progressively more rarefied as temperature decreases. 
Secondly, the resistances that comprise the CCN, as well as elemental ffuctuators that mod- 
ulate them, typically exhibit non-linear effects. At low temperatures, these non-linear effects 
are quite prominent even for small applied fields. As we shall see, these features introduce 
a new set of considerations into the question of noise correlations. 

In this paper we report on a study of noise statistics of crystalline indium-oxide films 
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{In20z-x) with their disorder tuned to make them strongly localized. Previous studies of 
1// noise using /„ 2 3 _ I fflms has focused on the raetalUc regime and employed macroscopic 
samples [8]. The emphasis in the current study is on the results obtained with samples that 
exhibit prominent mesoscopic effects. The other aspect that distinguishes this work from the 
previous study (where room temperature measurements were employed) is that the noise 
experiments were carried out at liquid helium temperatures. The original motivation for 
these experiments was an attempt to find the spatial extent of the correlation length in 
the glassy phase of In20^ x films by monitoring the noise characteristics of such films, in 
particular, the second-spectrum [9|, as function of their size. The rationale being that below 
a certain sample size the correlated nature of the electron glass should manifest itself in a 
respective correlation of the conductance fluctuations. That should presumably occur once 
the spatial scale associated with the glassy effects exceed the sample size. Samples with 
lateral dimensions ranging from 1mm down to 2/xm were studied for this purpose. A colored 
second spectrum, indicative of correlations, was indeed found for samples with sizes at the 
lower part of this range. However, the phenomenology encountered in these measurements 
was not in line with the expectation based on the above scenario. In particular, even in the 
small size regime, the occurrence of colored second spectrum was not consistent; samples 
with nearly identical parameters gave conflicting results. This led us to suspect that the 
source of correlations is not related to the underlying glass. Further experiments revealed 
that that the degree of coloration in the second spectrum depends in a non trivial way on 
the bias V used in the measurements. We propose a heuristic picture that qualitatively 
accounts for these observations. This model is generic to transport in disordered systems 
where transport is inherently inhomogeneous, and does not depend on the system being 
glassy. 

II. SAMPLES PREPARATION AND MEASUREMENTS TECHNIQUES 

The ln 2 0^ x films used here were e-gun evaporated on either a 110/im thick microscope- 
cover glass, or on a Si02 insulating layer (0.5/im thick) thermally grown on a Si wafer. The 
latter was boron doped and had resistivity 2 • 10~ 3 Qcm, deep in the degenerate regime. It 
thus could be used as the gate electrode for a low-temperature measurement with the sample 
configured as a MOSFET device in which a thin film of In^Oz-x served as the active layer. 



3 



Lateral dimensions of the samples were controlled using a stainless steel mask (for samples 
larger than 0.5mm), or optical lithography (for sizes in the range 30 — 200/im) and e-beam 
lithography for samples smaller than 30/xm. Samples used in this study had length (L) and 
width {W) that ranged from 2mm down to 2/xm and typical thickness d — 55 ± 5A The 
'source' and 'drain' contacts were made from thermally evaporated ~ 500A thick gold films. 
Fuller details of sample preparation and characterization are given elsewhere jiol ]. 

Conductance measurements were performed using two terminal ac technique, employing 
ITHACO-1211 current preamplifier and PAR-124A lock-in amplifier. In the MOSFET-like 
samples, gate voltage sweeps were affected by charging a 10/iF capacitor with a constant 
current source (Keithley 220). All the measurements reported here were performed with the 
samples immersed in liquid helium at T = 4.11K held by a 100 liters storage-dewar. This 
allowed long term measurements without disturbing the samples as well as a convenient way 
to maintain a stable bath temperature. These requirements are of particular importance 
for s t udl e S of g .a Ssy system where sample hrstory may influence time, dependent _ 
ments, as was demonstrated in previous studies using such samples [10]. To cater for this, 
samples were allowed to equilibrate for at least 15 hours prior to any conductance versus 
time measurement. No change in the nature of the noise of a given sample was found when 
measurements were repeated a week after the initial cool-down. 

A typical hopping length at these temperatures, and for the range of resistances used in 
this work is ~ 200A jljj, which makes such samples effectively two-dimensional (2D). This 
fact will be used in the proposed model described in section IIIC. 

Noise measurements employed a two terminal technique. These were performed by biasing 
the sample with a constant (dc) voltage source (home-made rechargeable battery-stack) while 
measuring the resulting current fluctuations by the ac voltage drop V ac detected across a 
series resistor. This V ac was monitored by a Dynamic Signal Analyzer (HP35670A) buffered 
by EG&G5113 low-noise pre-amplifier. These data were then used to calculate the first and 
second power-spectrum. The latter were implemented by the method suggested by Restle 
et al |l2|. As will be demonstrated below, the noise magnitude of the samples studied here 
is quite large. This is due to three reasons; (1) The samples reported here are strongly 
disordered and exhibit large noise even when of macroscopic size j^. (2) The samples 
are physically small, c.f., figure 2. (3) At the measured temperature, the samples are rather 
deep into the hopping regime, which means that the effective volume for the current carrying 
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process in these samples may be considerably smaller than their geometrical volume. These 
factors combine to give a large noise magnitude, which made it unnecessary to use a more 
comprehensive measurement configuration such as a 5-probe technique (that is awkward 
to employ for insulating samples being inherently inhomogeneous) . We did however check 
that the contacts do not contribute to the noise by measuring few samples using a 4-probe 
technique. 

The stability of the temperature bath afforded by the 100 liter storage dewar was suffi- 
ciently accurate to neglect time-dependent temperature effects on the samples conductance 
relative to the inherent 1/f noise. This was ascertained by monitoring the resistance fluc- 
tuations of a Ge thermometer attached to the sample stage as in Fig. 5 of Vaknin et al 

3- 

III. RESULTS AND DISCUSSION 
A. Preliminary measurements 

A conductance time-trace and associated power-spectrum for the smallest sample in the 
series studied are shown in Fig.l. Apart from its large magnitude (presumably due to the 
small sample size), the noise has similar characteristics as these of metallic samples of this 
material jq, in particular the overall shape of the power-spectrum retains its 1/ f n character 
although some deviations from a power-law may be observed; a best-fit to the data yields 
7] ~ 1.1 but the low / region seems to follow a faster dependence. This is not surprising as 
such samples contain a rather small number of fluctuating resistances as will be shown later. 
The dependence of the noise magnitude on sample volume is shown in Fig. 2 for several of the 
studied samples. To facilitate comparison, the figure includes only samples that are all of 
the same thickness, and have similar resistances. Note that the noise magnitude is inversely 
proportional to volume down to sample size of 2x2/im. This result is consistent with what 
is expected of ensemble-averaged wn-correlated fluctuators. It might therefore suggest that, 
if there are correlations in the noise in this system, they occur on scales smaller than 2//m. 
Surprisingly then, when analyzing preliminary data for samples in this series we encountered 
several traces that yielded a frequency-dependent second-spectrum. Such occurrences were 
not encountered in samples with sizes > 100/xm, and they were still rare in the range of 
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sizes 30 — 50 fim. Colored second- spectrum could be observed in samples as large as ~ 40/im 
but in unpredictable way; Two samples with nearly identical parameters (resistance, size), 
and measured under similar conditions, gave conflicting results; one exhibited frequency- 
dependent second spectrum, while the other had a Gaussian spectrum. An example of 
such a 'conflicting couple' is shown in Fig.3 for 30x40/im samples. This lack of systematic 
behavior persisted down to our smallest sample sizes, although the frequency of a colored 
second-spectrum appearances seem to grow with decreasing size. 

In contrast with the inconsistent appearance of noise correlations, all samples in this 
series systematically show the same glassy features as macroscopic samples with similar 
parameters. For example, Fig.4 shows the 'memory' cusp, which is the earmark of the 
electron-glass ^(|, for the 30xA0fim sample shown in Fig.3. Note that this actually is the 
sample that exhibits Gaussian noise. It is therefore unlikely that the /-dependent second 
spectrum can be related to the correlations due to the interactions that are associated with 
electron-glass dynamics. This conclusion will be re-enforced by a more elaborate analysis 
and the ensuing discussion below. 

In addition to the glassy cusp, the data in Fig.4 show mesoscopic conductance fluctuations 
(CF). These are the 'fingerprints' of the underlying CCN, and reflect the process by which 
some critical resistors in the CCN are replaced by other critical resistors as the chemical 
potential is varied Q|. The relative magnitude of the CF is a function of C, the correlation 
length of the percolation network. This is based on two assumptions: First, the basic 
conductance swing AG associated with a critical resistor is of the order of its conductance 
G. This is a characteristic feature of the strongly localized regime, (which for the present 
case of 2D samples means sheet resistance Rq that fulfills Rq 3> -§) Secondly, the 

relative fluctuation amplitude AG/G for the entire sample is essentially determined by the 



square root of the number of critical resistors in the sample N ~ therefore, ~ y j^- 
This relation will be used in this paper to estimate C as well as it dependence on applied 
voltage. 

B. Dependence on applied field 

Low temperature transport measurements on hopping systems are very sensitive to the 
value of the voltage V used. It is notoriously difficult to achieve linear response conditions in 
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such cases, especially for small samples at low temperatures, and deviations from Ohm's law 
are hard to avoid. This is illustrated in Fig. 5 for the two samples that will be discussed in 
this subsection. Usually, the effect of a not-small-enough-voltage on the measurement is in 
the same direction of raising the sample temperature. For example, the resistance decreases 
with V ( Fig. 5), and so does the noise magnitude (figure 7c), in both cases, the effect is 
monotonous with V. It turns out that the degree of non-Gaussianity of the type considered 
here, behaves in a qualitatively different way. In particular, it is non-monotonous with V, 
and it tends to vanish at both, high bias and at low bias thus peaking at an intermediate 
value of bias. 

As a quantitative measure of the degree of non-Gaussianity we take the area under the 
curve of the second spectrum as depicted in Fig.3, which we label as IntiyS?). Using this 
scheme, we plot the dependence of this quantity on V for one of our smallest sample that 
exhibits non-Gaussianity (Fig. 6). The first set of data (solid circles) were taken initially 
without regard to the order of changing V just to test the effect produced. When it was 
realized that the correlations may disappear at high V, the field was taken to a much larger 
value (V = 46 mV) to better define the asymptotic behavior of Int(S2)- A later attempt to 
add more points to the curve failed to reproduce the position of the peak in Int(S2) versus 
V obtained in the first series. Rather, the curve seems to have shifted towards a lower bias. 
The average sample resistance r and its dependence on voltage r(V) were not affected by 
the high V exposure. We shall return to this 'mesoscopic' effect following the discussion in 
the next section. 

Figure 7 shows a more elaborate study of another 2x2/j,m sample with similar noise 
characteristics except that now care was taken not to subject the sample to an excessively 
large V. The same pattern in terms of Int(S2) versus V emerged (upper graph in Fig. 7), and 
this time the curve was fairly well reproduced by a second series of measurements. Along 
with the second spectrum analysis, the figure shows the respective dependence of the power- 
spectrum parameters r\ and amplitude. Neither shows the non-monotonic V dependence 
exhibited by Int{S2)- It is natural that above a certain field, both Int{S2) and the noise 
magnitude decrease with V; a large field, like temperature, decreases the hopping-length and 
C, which in turn means larger number of fluctuators to average over. The tale-telling result 
is that, below certain field, Int(S 2 ) diminishes when V decreases in a way that suggests a 
much smaller effect as V — > 0. In other words, the noise appears to be correlated because 
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the measurement is not in the linear response regime. 

A plausible scenario that leads to correlations between individual fluctuators is based on 
two indents; A) the eurrent ta the Sy8 te m is earried by a percolation-— BflQ, 
and B) the (average) frequency Ui of a fluctuator i is, among other things, a function of 
the local voltage V^. The latter is true for a generic two-level-state whether of 'atomic' or 
'electronic' nature; The frequency of the two-level-state usually depends exponentially on the 
local voltage Q], so this is a sensitive source of inter-modulation once the local V^'s are re- 
distributed by the DCR effect. Ingredient A is an inherent feature of variable-range-hopping 
systems, and it is probably a common feature in other disordered conductors as well. The 
system may then be viewed as a random- resistor- network where the active fluctuators are 
part of it (or situated nearby such that can modulate a resistance that is in the current path). 
When a fluctuator j in the network, changes its state, the local voltage on fluctuator i will 
change too due to the continuity of the current carrying network. If the resulting voltage 
change 5Vi is not much smaller than fc^T, the switch of j will result is a change Sui of 
fluctuator i frequency. Such an effect is the basic building-block of a hierarchical correlation 
chain where the slower fluctuator modulates a faster one. Naturally, this mechanism for 
non-Gaussianity must vanish with the applied field as indeed is observed. It should also 
be negligible when the system size is much bigger than the range of the proposed inter- 
modulation effect. These considerations will be now dealt with in a more formal way, and 
the results will be compared with the experiments. 



C. Theoretical considerations and comparison with experiments 

We consider a 2D hopping system and use the standard percolation scheme [n| . Focusing 
attention on a given critical resistor j which is affected by some fluctuator. For simplicity 
we assume that each resistor is coupled only to a single fluctuator and will denote this 
fluctuator by the same index j . Two different candidates for the role of the basic fluctuators 
may be considered. The first is an 'atomic' two level system as in amorphous materials 
(see, e.g., Mott and Davis [17]). The second one is an aggregate composed of localized sites 
(not necessarily part of the percolation cluster), and having two metastable configurations 
characterized by different distribution of the electrons over the sites. The simplest object 
of such a sort is a pair of hopping sites occupied by a single electron considered as a source 
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of 1/f noise in jig ]. Although quantitatively the effects of atomic and electronic two level 
systems are expected to be different and depend on different parameters of the material, 
in both cases the fluctuators are expected to be sensitive to the local electric fields (for 
the atomic two level systems the corresponding coupling is related to a presence of electric 
dipolar moment of the atomic two level system). The model considerations given below 
are applicable for either mechanism, although the microscopic equations for the coupling 
coefficients are naturally different. 

The fluctuation of the corresponding voltage swing, 5Vj, will inevitably lead to a fluctua- 
tion of the voltage across all other resistors. This is the DCR effect considered by Seidler et 
al P and shown to lead to a colored second-spectrum of magnitude S^dor- Next we show 
that in the hopping system the voltage re-distribution in the CCN results in a more elaborate 
coupling mechanism between different fluctuators. In particular, this results in a dependence 
of the switching rate of the fluctuator i on the state of the fluctuator j. This coupling gives 
an additional, non-linear, contribution to the frequency dependent second spectrum (S2, C orr)- 
At large bias both S2,dcr and S2, C orr are suppressed with V due to an increase of the number 
of effective fluctuators N. We will show that S2 }C orr initially increases with V and peaks at 
some intermediate V consistent with the experimentally observed behavior (figures 6 and 
7). 

Denoting by SVji the voltage swing on resistor-i affected by the corresponding variation 
of resistor-j, and assume that the main ensuing effect is a variation of Aj, the difference in 
energy between the two states of the fluctuator; 

6(Ai) = B t 5V jt (1) 

where Bi is a coupling coefficient. In terms of the correlation length C, and for ^> C 
(Rij is the distance between the critical resistors % and j) SVji can be estimated as: 

sva - (2) 

(jj no- 
where in the regime eVj « ksT one has \SGj/Gj\ = \5ej/kBT\, and Sej is related to the 
modulation of the resistor activation energy by the current re-distribution. 

The variation of Aj leads to a variation of the fluctuator dwell times. The latter can be 
written as j3] 

r 1 = ^ + rJ = (l + = r-\l - n)^ (3) 



where r + and r_ are the dwell times for the upper and lower level, respectively, and n is the 
occupation at the upper-level state. Note that r + corresponds to a transition accompanied 
by the emission of the phonon with a frequency uj = A/h. Thus, the change in r + due to 
the fluctuation of A is 

5t + ~ -r + ^(SA) (4) 

where we have taken into account that T+ 1 oc A a . Note that a ~ 3 for either a fluctuator 
of electronic nature in the limit of small A, Q| , as well as for a fluctuator of atomic nature 

Combining eqs. I3l4l we obtain for the fluctuation of the relaxation time of z-th fluctuator: 

*r< = -n(^~ 5A,, (5) 



.A, k B T, 

Note that for an ideal 1/ / first spectrum (which, strictly speaking, can be realized only for 
L — > oo), the fluctuations of Tj would not lead to a variation of S\. However for a small 
sample size (where deviations from 1/f spectrum may be observable as, e.g., Fig.l), the 
first spectrum may be significantly affected by fluctuations of Tj. In the limit eVj << k B T 
(where £ = C ) the effect is proportional to Vj and thus to the total bias V (Vj ~ 

It can be shown (see Appendix 1) the resulting contribution to £2 can be estimated as 

S 2 ,corr OC 0?Bi8Vij OC V (6) 

By comparison, the contribution of the DCR is given as 

S 2 ,dcr oc Tjij (7) 

Here 



5V 2 



which at eVJ << k B T does not depend on V while 7^ means an average magnitude of 7^; 
the prefactors in eqs. 16171 differ from one another only by a numerical coefficient of the order 
of unity. 

Eqs l6l7l then suggest that S 2 , C orr > S 2t DCR if 

eVi > k B TJ^- (8) 
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In other words, our mechanism dominates over the DCR even when eVi <C ksT provided 

lij < o?B % (9) 

Recall that Bi describes a relative effect of the fluctuation of the resistor potential on the 
state of the nearby fluctuator while 7^ describes the relative fluctuation of the voltage on the 
resistor i due to a fluctuator modulating resistor j. In general, Rij > C and this ratio is small. 
At the same time the effect of the resistor voltage on the fluctuator can be large enough. 
It holds in particular for the fluctuators of the electronic origin if their size is comparable 
to the inter-site distance within the hopping resistor provided that it is situated at the 
distances less or comparable to r^. It also holds for the structural fluctuators provided they 
are situated close enough to the site with a lower energy. In these cases Bi ~ 1 therefore 
condition El may be obeyed. 

Now let us discuss the regime Vi > ksT when the hopping conductivity is strongly non- 
linear. It can be shown that the fluctuation of the conductance of resistor i resulting from 
the fluctuation of its activation energy for eVi > k^T is still given by dGi ~ — G^de^/fc^T. 
However the relative fluctuation of the total conductance of the sample depends on the 
correlation length C: 

o t \ — (dG^dG)^ _ C A (dGi, dGi) nn -, 
Sl{uj)= G 2 -tf^> Gl (10) 

where L is a linear size of the sample while £, is the correlation length of the percolation 



21] 



cluster which in the nonlinear regime can be estimated as 

where E is an average electric field within the sample, Cq = C(V — > 0), and v is the 
percolation theory index (for 2D v ~ 4/3). Thus £ oc V~ u ^ 1+ ^ = V~^ 7 . Correspondingly, 
one has 

Sl (tu)oc^ocV- 8 / 7 (12) 

(we have taken into account that the result of a summation over the effective resistors is 
proportional to the number of these resistors). 

This ensemble- averaging effect is expected to be even stronger on the second spectrum 
since it is a convolution of the two first spectra. So one expects S2 oc N~ 2 where N = (L/C) 2 
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is the number of fTuctuators, then: 

S2 ' corr K N 2 (V) ^ 
This is strictly obeyed for eV < k B T, at higher voltages there may be contributions of other 
non- linear mechanism not considered here. Note that Eq.13 contains two factors. The 
numerator (V) depends on the coupling coefficient given in Eqs. 1-2 (see Appendix) and is 
indeed defined for the case eV < k B T. For larger bias this coupling coefficient is suppressed 
due to the non-linearity of the medium and corresponds to a sub-linear behavior; we do not 
consider this effect in detail. The denominator depends on the statistical average discussed 
above Eq.13, and it remains the same for eV > k B T while its dependence on V follows from 
the considerations given in Eq.ll. 

To compare, the respective contribution of the DCR mechanism has the following depen- 
dence on V : 

S2,DCR OC (14) 

Let us now see how these expectations compare with our experiments. To find the 
qualitative dependence on V of the second-spectrum amplitude we need to know the function 
N 2 (V), namely how the number of fluctuators varies with V over the range relevant for the 
experiment. This may be estimated theoretically using similar considerations as those that 
led to eq^J above. The dependence of the first spectrum amplitude on V is in rough 
agreement with this equation (c.f., Fig. 7c). However, the data for the second spectrum 
amplitude were taken over more extensive range of V, exceeding the limits of validity of the 
power-law relation expected by eq!12l It is therefore necessary to get an estimate for N 2 (V) 
from experiments. That can be done through the use of data for the conductance fluctuation 
versus V such as the results shown in Fig.8a. To construct an empirical N 2 (V), one then 



uses the relation ^p(V) ~ y = N(V) 2 discussed in section 2. This procedure yields 
the N(V) depicted in Fig. 8b, which empirically, fits rather well an exponential dependence; 



N(V) oc exp(vF ) 22j. Using this form in eqs ll3ll4l one gets the qualitative dependence on 
V for the two mechanisms for the second-spectrum considered here. This is schematically 
illustrated in Fig.9. The overall shape in this plot is in fair agreement with the experimental 
curves (c.f., Fig. 6 and 7) although it seems that the small bias regime would fit better a 
faster than linear with V relation. To estimate the value of the voltage Vi (that at the peak 
of Int(S 2 ) versus V should probably be of the order of k B T, c.f., Fig.9), one needs to know 
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the number of critical resistors in the sample. As noted above, this can be done based on 
the relative magnitude of the CF, namely, the reproducible fluctuations in G(V g ) generated 
by sweeping the gate voltage V g . The number of critical resistors along the sample in Fig. 6 
can be be estimated using the data in Fig.8b. Note that the applied voltage at the peak of 
Int^S?) is V^, ~ bmV and V p ~ 10m V for the data in squares and circles respectively (c.f., 
Fig. 6). The value N for V ~ lOmV can be read from Fig.8b as N ~ 260, which then means 
that the value of Vi at the peak is ~ ~ 0.6mV. Similarly, The value N for V ~ 5mV 

is iV »s 91, giving Vi ~ ~§Q m V ~ 0.52mV. These values compare favorably with the sample 
temperature T = 4.1 IK. 

We were not able to measure noise in the 2/xm samples using bias that is strictly in the 
linear response regime (eV, <C kgT). In fact, deviations from linear response in the range of 
bias used here are reflected in the sample conductance itself (Fig. 5). The smallest bias used 
for the sample in Fig. 6 was V = lAmV which, in terms of V,, is the equivalent of 3K (using 
the respective N ~ 29 from Fig. 8. This bias is not much smaller than the bath temperature. 
Nevertheless, the low bias regime we did manage to use is low enough to expose the peak in 
Int{S2) vs. V consistent with the proposed mechanism. 

Note that the relevance of both correlated-noise scenarios discussed here hinges on spe- 
cific assumptions. The DCR assumes that the conductance swings associated with slow 
fluctuators are potent enough to give a significant contribution. The non-linear mechanism 
we offered requires the existence of 'soft' fluctuators that, in addition, are coupled effectively 
to critical resistors. In either case the 'master fluctuator(s)' should operate on the frequency 
window that is relevant for the experimental scales. In the regime of mesoscopic samples 
one may expect to find fluctuators realizations such that some (or all) of these conditions 
are not realized in which case the non-Gaussianity will be weak or absent. Moreover, ap- 
plying a voltage may displace a key fluctuator out of its 'commanding' position thereby 
removing the origin of correlations. This is certainly a concern in the strongly non-linear 
regime. Applying a large field will inevitably modify the current-carrying network. The 
modification may be reversible, in which case we expect that data such as in figures 6 and 
7 will reproduce themselves under V cycling. However, when the applied V is sufficiently 
large, it is quite likely that a different percolation network will be precipitated, just as a 
thermal cycling involving high temperatures would cause j23j. This may lead to noise data 
of a different statistical nature. The change may well be subtle; turning on or off the cou- 
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pling of certain fluctuators to the CCN is all that is needed. The average disorder may not 
change in the V - recycling process, and the CF pattern may be only slightly affected. We 
believe that these considerations account qualitatively for the non-systematic occurrences of 
non-Gaussian noise in our samples, as well as for the effect of V cycling discussed in section 
II. 

It should be mentioned that there are other non-linear mechanisms, not considered here, 
that may contribute to correlations between remote fluctuators, especially at the high V 
regime. For example, the change of the local voltages could modify the values of the critical 
resistors in the CCN j^j], and may give rise to new fluctuators. Another interesting scenario 
is a heterodyne effect; namely, frequency mixing of the 'master' frequency with the 'local' one 
due to the non-linearity of the critical resistors. All such mechanisms, as well as the DCR, 
should be seriously considered whenever a noise with colored second-spectrum is encountered 
in a conducting system. Note that a generic feature of these mechanisms is their long range 
1/R nature (see eq|2]) making them more effective than most 'direct' interactions. As was 
demonstrated in this work, these effects may bring about correlations between fluctuators 
even in samples that are considerably larger than the scale relevant for the interactions 
associated with the electron-glass j^ . 

Finally, it should be mentioned that while reducing the applied V below the value where 
Int(S2) peaks diminishes the non- Gaussian effect, the glassy effects if anything, become 
more prominent HQ- Clearly then, glassmess and correlated -noise (when exists) are not 
necessarily related. 

In summary, we have described a set of conductance-noise experiments on disordered 
films of In 2 Os- x in their glassy phase. The emphasis in this study was on the degree of 
noise correlation as function of system size. Noise correlation was measured by the second 
spectrum of conductance data taken at liquid helium temperatures. Our main finding is 
that, down to sample size of 2/im the noise has the usual 1/f power-spectrum. Hopping 
samples with this size contain small number of critical resistors as indicated by the prominent 
conductance-fluctuations they exhibit. Such samples still show essentially all the electron- 
glass features as macroscopic samples j^. Given the way non-Gaussianity decreases below 
a certain bias in these samples it seems unlikely that the correlations observed at finite bias 
are due to glassiness. 

When the voltage used in the conductance measurements was not small enough, non- 
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Gaussian noise was observed in several of the samples, including samples as large as 40/im. 
We have demonstrated that the degree of non-Gaussianity is a non-trivial function of the 
bias. A model that purports to account for these findings was offered, and its consequences 
are found to be in qualitative agreement with our experiments. 

A lesson that may be taken from our study is that correlations in conductance noise may 
arise from a non-linear mechanism, and this may be of particular relevance to disordered 
conductors measured at low temperature. Such effects need be better understood and care- 
fully examined before a non-Gaussian noise is associated with correlations due to, e.g., glass. 
As a point of principle, one indeed expects that, on sufficiently small scales, a glassy system 
may show a correlated noise. What should perhaps be stressed is that the converse is not 
necessarily true. 

We acknowledge illuminating discussion with Clare Yu on her computer simulation results 
of the DCR model. One of us (V.I.K.) acknowledges support of the Lady Davis Foundation. 
This research has been supported by the Binational US-Israel Science Foundation and by 
The Israeli Academy for Sciences and Humanities. 



IV. APPENDIX 



Here we consider the coupling between the ffuctuators in more detail, and estimate the 
effect of this coupling on the second spectrum. For simplicity we assume that the (initially 
un-correlated) 'fast' ffuctuators are coupled to 'slow' ones, and the coupling is characterized 
by the occupation numbers hj. By 'fast' and 'slow' ffuctuators we refer to specific two- 
level-systems that contribute in the measured noise, (at the high end of the ffuctuators vs. 
frequency distribution and at the lower end respectively). Explicitly we consider only the 
modulation effect of 'slow' ffuctuators on 'fast' one. The complementary process (i.e., 'fast' 
affecting 'slow') while possible is much more involved and is not treated here. 

The 'slow' ffuctuators affect the value of Aj of the 'fast' ffuctuators. A natural result of 
this modification is change of the ffuctuators occupation numbers, rij. 

Making use of Eq|TJ and using n = (exp(A/fc B T) + 1) _1 one concludes that in a presence 
of 'slow' ffuctuator j one has 

n,(l - n0|, = m{l ~ nMl - tanh ^/^ T ) ' ^.(t)^) (15) 
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Another source of the fluctuations is related to fluctuations of the relaxation time T;. One 
readily obtains: 

6 ( tTT vr l = s ^ 7TTT%i (16) 

Thus, the contribution of i — th resistor to the first spectrum is: 

5 ((5m, (Jni) w ) = Aijhf (17) 

a /r r \ i tanh(Aj/2/;;BT) _ _ Tr 
Aij = -(dTii,dni) w | jrp-^ rij^BidVij- 

- t - (i + m W "- (1 -"- ) U'^J w " 

As a result, the first spectrum appears to be dependent on the occupation number of the 
'slow' fluctuators. 

Now let us estimate the effect of the correlations on the second spectrum S^u^). By 
definition for S 2 we have 



S2(w 1 ,w 2 ) = G- 4 / dre iu)T / dt 
Jo Jo 

»t+t pt'+t 

±/\ it'uji 



dt' J dT'5G(t' + T')5G(t')e l 

ft+T+to ft"+t 

/ dt" / dr"5G(t" + r")5G(t")e lT "^ (18) 

Jt+T Jt" 



where in our case 



6G = Y^9iSrn (19) 



Thus, the spectrum of the fluctuations is related to the temporal behavior of Srii. Here are 
the coefficients describing a contribution of z-th fluctuator to the conductance fluctuations. 
Let us first consider a case of statistically independent fluctuators. In this case 

5G(t')5G(t' + r') = J29i( s ni(t')Sn(t' + O) (20) 

i 

The DCR mechanism results from the modulation of the coefficient gi in Eq^Jby another 
('slow') fluctuator j. The corresponding contribution to S 2 can be estimated as 



S2,dcr oc jjj(8ni, Srij)^ (dnj, 5r 

Pk B Tf 2 (ui)% ^2(5nj, Sfij)^ (21) 
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In contrast, in our case we deal with real coupling between the fluctuators described 
by Eq^ which exists only at finite V. Using EqEI one obtains a contribution to S 2 
oc Afj(5hj,5hj). The most important term here is the one oc (a/Aj) 2 . Indeed, due 
to a presence of 'soft' fluctuators with Aj — > this term is divergent. The divergency has a 
clear cut-off at Aj ~ <5Aj = B^Vij since in our calculations we assumed the fluctuations of Aj 
to be small (note that here we assume eVj < ksT). Summing over the different fluctuators 
gives 



where we have assumed that the distribution -P(A) = (P) is flat. Thus here the summation 
over different A is controlled by the lower limit of 'soft' fluctuators, and we obtain 



3 

Assuming that the 'slow' fluctuators have exponentially broad scatter of relaxation times, 
the contribution to the second spectrum is S 2jCorr oc uj 2 x (at the frequency scale t$ 1 > co 2 > 



Note that the non-linear effect we propose will dominate over the DCR mechanism even 
for moderate bias eVj < ksT. The reason is that among different fluctuators there exist those 
with small Aj. For such fluctuators the relative modulation of the relaxation rate due to 
fluctuations of Aj as it is seen from Eq. Elappears to be oc 5Aj/Aj. Since the corresponding 
contribution to S 2 is proportional to (5Aj/Aj) 2 it leads to a significant enhancement of a 
role of 'soft' fluctuators with small Aj as is indicated by eq. 1221 
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V. FIGURES CAPTIONS 



FIG. 1. Conductance noise power-spectrum measured under bias voltage V = 3mV. The 
lower frequency part of the spectrum represents the averaging over 16 Fourier transformed 
G(t) runs. The higher frequency part of the power-spectrum was averaged over 1024 power- 
spectrum curves measured on the resistor R = 1.01MJ1, connected in series with the sample, 
employing HP35670A. The dashed line depicts 1/ f im dependence (a best fit to the data). 
Inset: typical G(t) series used for calculation of the lower frequency part of the power- 
spectrum. Note the relatively low signal to noise ratio typical for mesoscopic samples. The 
sample: length=2/im, width=2/im, R\j = 5MQ. 

FIG. 2. The magnitude of the noise power-spectrum measured at / = lOmHz as a 
function of sample area A (= LxW). Each datum point represents the averaging over 
several samples with same area and values of the resistance ranging between Rq = 2MQ 
and R\j = 10Mf2. The dashed line depicts the Si oc A' 1 law. 

FIG. 3. Second-spectrum as a function of normalized frequency measured in two samples 
with the same lateral dimensions; L = 30/im, and W = 40/im, and resistances; R\j = 
4.8MJ1, and Rq = 5.3Mfi for samples in (a) and (b) respectively. The octaves (2 to 9) 
are labelled by the corresponding lowest frequency values (/x,) and represented by different 
symbols on the plot. The measurements for both samples were performed using a series 
resistor R = 100/cfi and applying V = lOOmV, and V = 150mV as a bias voltage for the 
sample in (a) and (b) respectively. The degree of non-Gaussianity (labeled in this work as 
Int(S 2 ), see text) is taken to be proportional to the area defined by the dashed and dotted 
lines, and the ordinate axis (plate a). 

FIG. 4. The 'memory' cusp (see reference 9) as seen in the measurements of conductance 
as a function of the gate voltage for the sample shown in Fig.3(b). Two successively 
measured traces (solid and open circles) show reproducible conductance fluctuation (CF). 
Gate voltage scan rate was 0.02V /sec, bias voltage V = 20mV. 

FIG. 5. Conductance as a function of the applied bias measured for two samples with 
the same lateral dimensions; L = 2/im, W = 2/im. Note the deviation from linear response 
even at the smallest bias used. 



19 



FIG. 6. The degree of non-Gaussianity Int(S 2 ) as function of the applied bias for two 
series of measurements: prior to (solid circles), and after (open squares) the application of 
high bias (see text). Dashed lines are guides for the eye. The measurements were performed 
using a series resistor R = 1.01MQ. The sample: L = 2/j.m, W = 2/zm, i?n = 5MQ. 

FIG. 7. Plate a: The degree of non-Gaussianity as function of the applied bias in two 
series of measurements: successively increasing the bias (open circles), and a later set, 
employing bias values within the same range as before with no particular order (solid 
circles). The value of rj (in the first-spectrum law, 1/ f v ), and the noise power per decade 
as function of the applied bias - (b) and (c), respectively. Dashed lines are guides for 
the eye. The measurements were performed using a series resistor R = 1.01MQ. Sample 
parameters: L = 2/im, W = 2/im, Rq = 9M£l. 

FIG. 8. Plate (a): Conductance versus gate- voltage scans taken with different values of 
bias V for the sample in Fig.6. Each trace shows reproducible conductance fluctuations with 
rms amplitude that decreases with V. The average value of the conductance increases with 
bias (c.f., Fig.5). Plate (b): The values of N for the series of traces shown in (a) estimated 
by the rms amplitude of the fluctuations (see text). 

FIG. 9. A schematic dependence of S 2 , C orr and S 2 ,dcr on the applied voltage (see text). 
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